Probing Dark Energy with Black Hole Binaries 
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The equation of state (EoS) of dark energy w remains elusive despite enormous experimental efforts to pin 
down its value and its time variation. Yet it is the single most important handle we have in our understanding of 
one of the most mysterious puzzle in nature, dark energy. This letter proposes a new method for measuring the 
EoS of dark energy by using the gravitational waves (GW) of black hole binaries. The method described here 
offers an alternative to the standard way of large scale surveys. 

It is well known that the mass of a black hole changes due to the accretion of dark energy but at an extremely 
slow rate. However, a binary of supermassive black holes (SBH) radiates gravitational waves with a power 
proportional to the masses of these accreting stars and thereby carries information on dark energy. These waves 
can propagate through the vastness of structure in the universe unimpeded. The orbital changes of the binary, 
induced by the energy loss from gravitational radiation, receive a large contribution from dark energy accretion. 
This contribution is directly proportional to ( 1 + w) and is dominant for SBH binaries with separation R > 1000 
parsec, thereby accelerating the merging process for vv > — 1 or ripping the stars apart for phantom dark energy 
with w < — 1. Such orbital changes, therefore w, can be detected with LIGO and LISA near merging time, or 
with X-ray and radio measurements of Chandra and VLBA experiments. 

PACS numbers: 98.80.Qc, 11.25.Wx 



I. INTRODUCTION 



One of our most crucial questions about nature at present 
is: what is dark energy? The fact that our universe is acceler- 
ating [1, 3] and that dark energy constitutes about 70 percent 
of the total energy density of the universe, are well established 
by now. Many theoretical models have been put forth which 
cast dark energy in the form of a cosmic fluid tH 111, It'] with 
time variations in its equation of state w{z)- Yet the simplest 
explanation for dark energy remains to be a pure cosmological 
constant (cc) A. 

The trouble we face in understanding dark energy does not 
stem from a shortage of dark energy models, with w (zo ) — — L 
that mimick at present the behavour of A and give rise to the 
observed acceleration of the universe. The puzzle rather lies 
on identifying which one of these possible candidates is the 
correct one. The best way to discriminate among the various 
possibilites and a pure cc, A, is to experimentally measure 
the time variations of the dark energy equation of state w(z). 
So far a popular parametrization for w{z) is the linear one 
w(z) = Wo + wiz+ ... with wi = for a pure cc A knowl- 
edge of w{z) is crucial for not only understanding the present 
accelerated expansion of the universe but also for making pre- 
dictions for its future evolution and destiny. If w{z) "C — 1 
then dark energy is a phantom |9] which leads the universe to 
a Big Rip in the future. If w = — 1 then we are probably ITlll 
facing an eternal DeSitter state 1 12] which at least at the clas- 
sical level implies constant temperature and entropy therefore 
a cosmic heat death \2l, 22]. Other forms of w{z) can also 
allow for a Big Crunch 1 10] or bounces ] 17]. At present we 
can not infer which destiny our universe will meet without a 
better knowldege of w{z)- 
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Major experimental efforts for pinning down the value and 
time-variations of w{z) are under way through large scale sur- 
veys from CMB ]1], large scale structure and SNla ob- 
servations. The endeavor of measuring, to confident preci- 
sion, such small time variations in w{z) has proven extremely 
difficult, partly due to the inherited errors in the experiments 
that are not instrumental but which originate from noise accu- 
mulated from the background and foreground effects through 
which the signal we receive has propagated. In order to mini- 
mize such errors associated with the propagation of the signal 
through the vastness of structure in the universe, we would 
like to propose in this letter a complimentary method for mea- 
suring w{z). This method uses compact and localized objects, 
such as Black Holes, for acquiring information about w{z), 
by exploiting the gravitational waves these objects emit when 
they are in binaries. The advantage of this method is twofold: 
first, gravitational waves propagate undisturbed through struc- 
ture; and second, we have existing experiments which are ei- 
ther already operational or will be in the near future, such as 
LIGO and LISA missions designed to detect these binaries 
gravitational waves, or Chandra and VLBA experiments de- 
signed for X-ray and radio measurements. 

The accretion of dark energy by Black Holes is reviewed 
in Sec. 2., including a review of the main parameters of 
binaries and gravitational waves, useful for our purposed. 
Sec.3presents the method we propose, along with an investi- 
gation, discussion and some illustrations, on how gravitational 
waves from SBH's binaries can be used for extrapolating the 
equation of state w{z) of dark energy. 



II. BLACK HOLE BINARIES 

The phenomenon of accretion of dark energy by black holes 
is now established. If dark energy is modeled as a background 
cosmic fluid, then the flux of energy accreted by the hole will 
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change the mass of the black hole M by a rate given by the 
equation ifllfll 



(2.1) 



where A ~ 5.6 is a numerical factor whose value depends on 
the matching done for the dark energy fluid velocity at the 
horizon of the black hole, poo is the energy density of dark 
energy far away from the hole, and w{z) is the equation of state 
for dark energy as a function of the redshift z- For the case of a 
quintessence field, Eqn l2.1l becomes M = ±47t(2M)^(j)^ where 
the (— ) sign corresponds to phantom fields[13]. Clearly by 
accreting dark energy black holes gain in mass for vv > — 1 
and lose mass for w < — 1 . The solution to Eqn l2.1| is 



M = 



Mo 



(2.2) 



where Mq is the initial mass of the black hole and T is its evo- 
lution time scale 
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47tAMoPoo(l +w) 



(2.3) 



In the case of phantom energy, T provides the time scale to 
the Big Rip, by which time all black holes must have com- 
pletely evaporated. 

If dark energy is a quintessence scalar field with w > —I, 
the authors of lfl4ll have shown that the growth of the mass of 
the black hole due to accretion of dark energy could provide 
a sufficient mechanism for converting primordial black holes 
(PBH) into supermassive black holes (SBH) within a reason- 
able time, a time comparable to the age of the universe. SBH's 
can be produced via dark matter accretion 1 16] and they can 
have a mass as large as where is a solar mass. 



A. SBH Binaries 

Most active galactic nuclei (AGN) are populated by SBH's 
therefore we can expect that SBH binaries are not that uncom- 
mon. A binary of two massive objects will emit graviational 
waves with angular frequency Q. 



G{mi +m2 
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(2.4) 



where R is the separation of the stars in the orbit, G is New- 
ton's constant, and {1111,1112) are the masses of the stars in the 
binary. The period of rotation is T = Another useful de- 
tection parameter is the amplitude = '^'"j^"'^ where r is the 
distance along the line of sight to the observer. 

The power Pqw contained in the gravitational waves results 
in a loss of the energy of the binary E — Mc^ ifTsll where 



Pgw — 
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R^ 



(2.5) 



and the effective gravitational mass of the binary M is de- 
fined by 



M = m\+m2 



1 m\m2 

2 R 



(2.6) 



As a result the orbits get smaller until finally the merging 
occurs. LIGO and LISA experiments are designed to look for 
such gravitational waves, while VLBA and Chandra experi- 
ments can observe the binaries via radio and X-ray measure- 
ments. The optimal LIGO band for detection is in frequencies 
around i2=150i/z in and LISA will be able to detect waves 
with frequencies as low as IQ^^Hz- 



III. INFORMATION ON DARK ENERGY FROM BLACK 
HOLE BINARIES 

In order to illustrate the main idea of the method pro- 
posed here, let us for simplicity take equal mass binaries, 
m\ = 1112 = in of supermassive black holes that are in the back- 
ground of the unknown dark energy fluid. In this case the ef- 
fective mass from Eqn. |2.6| is: M — 2m — and the expres- 
sion Eqn |2.5| for the power of the emitted gravitational waves 
(GW), becomes 



Pg 
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w 



2m^ 



R^ 



(3.1) 



Since the binary is immersed in the 'bath' of dark energy 
with energy density pA, the mass of each hole in the binary 
will change according to Eqns |2. 1112.21 due to the accretion of 
dark energy 



rii — m"^47tApA(l + w) = 2m^L 



(3.2) 



Here w{z) denotes the equation of state of dark energy at red- 
shift z, and mo the mass of the star at initial time f = 0. The 
parameter L is defined such that: 2Loto(^) = t^^, with the 



evolution time x given in Eqn. 12.31 If the mass of black holes 
in the binary changes then the power of the emitted gravita- 
tional waves Pew has to change accordingly, due to the accre- 
tion of dark energy. But the energy losses from Pqw prior to 
the inclusion of dark energy accretion and after, result in quite 
different orbital change and in some cases in a highly differ- 
ent merging timescale, (see Fig.l and Fig.2). These orbital 
changes due to dark energy depend sensitively on Pa,w(z),oto 
We can estimate the dependence of /? as a function of w and 
its time evolution, /?[r,vv] by solving a differential equation 
which is derived by matching the mass changes [dMc^/dt] 
which now are due to both: dark energy and GW's, with the 
energy losses Pqw through gravitational waves. 



Pgw — 



d{Mc^) _ -32 
dt ^ 5 



2m^ 



R^ 



(3.3) 



The effect of dark energy accretion on the mass of stars 
Eqn. I3.2l is now included in the expression for power Eqr j3.1l 
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and the matching in Eqn. I3.3l resuhs in the following differen- 
tial equation 
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2mQ 



dt 5 
-ALR^niQ 
[1 - ILtmo] 



[1 — 2Lfmo]- 



%LR^ 



(3.4) 



This equation, along with its solution below, Eqn. 13.71 are 
the most important results of this investigation. Notice that the 
first term for m — mo corresponds to the well known general 
relativity effect of gravitational waves emission, the Hulse- 
Taylor effect. But the other two terms are new and completely 
due to the dark energy accretion by black holes. 

Clearly the orbital separation /?[r,w] depends sensitively 
on the dark energy equation of state w(z). We can thus use 
observations of the orbital changes R with time in order to 
extract information about dark energy's EoS w{z) by using 
the already existing and planned GW's experiments, LIGO 
and LISA, or directly via X — my measurement missions like 
Chandra and radio measurements with VLBA ll20ll . 

The changes in the frequency Ci and therefore the orbital 
period T — ^ can be derived from Eqn l3.3l in the same man- 
ner, which yield 



For this reason, we focus on the large mass limit below and 
show that the dark energy correction terms to the orbit, Eqn. 
13.71 are significant and even dominant for SBH's binaries with 
Rq ~ \QQQpc,a~ 10^ or larger 

All the information about dark energy is contained in L = 
2^ = 27tpA(l + w) which is positive for all dark energy mod- 
els with w > — 1 . For the case of a pure cosmological constant 
(cc), L = 0, and the corrections due to dark energy disappear 
entirely. However for phantom dark energy where w < — 1 , L 
is negative. The correction term due to dark energy in Eqn l3.6l 
in the case of phantom dark energy, thus has the opposite sign 
to the third (Hulse-Taylor) term that originates from power 
losses due to emission of gravitational waves. The sensitivity 
of the binary to the dark energy EoS w{z) is now clear: orbital 
changes of a binary immersed in phantom energy L < in- 
crease the separation and, are quite different from the orbital 
changes in the case of dark energy with L > which accel- 
erate the merging, or the case L = where no changes of the 
orbit are induced from dark energy. As shown below, this re- 
sult is used to deduce whether dark energy is a cc, a phantom 
or whether w lies above the cc boundary w = — 1 . 

In order to get a family of examples, let us parametrize the 
masses of SBH's and the orbital separation by the following 
relations 



Rq 

R[t,w] 



(3.5) 



where Rq is the initial separation and £2o initial angular fre- 
quency at f = 0. If the orbits get smaller with time, ^[f , w] <C 
Rq, then according to Eqn. 3. 5. the frequency £1 might increase 
sufficiently such as to fall within the detection limits of LIGO 
or LISA. From solutions to Eqn l3.4l we can get an estimate for 
the effect of dark energy on parameters of the SBH's binary. 

To a first order approximation the solutions of Eqn |3.4| is 



R[tM = 
32 16 



[rI + 8/?^Log(l - ILtmo) - 32LRlt 

2 

,1/4 



(1 -2Lfmo)2 



In general T is quite large. For example, for a solar mass 
M*,T ~ lO^^yrs. Therefore to a good approximation, we can 
re-write Eq l3.6l including the appropriate units, as 



(3.9) 



where the parameter a quantifies how heavy the SBH is rel- 
ative to the solar mass M*, and 



Rq = 2mop(- 



:2PfllO^ 



(3.10) 



where the parameter p ^ 2 reinforces the requirement that 
the orbital separation better be larger than the Schwarzchild 
radius 2mQ of each star at initial time f = 0, long before merg- 
ing. 

Since the evaporation time scale x for a solar mass black 



(3.6) hole (fl = 1) is about T — 10 s = 10 yrs, then the Ufetime x 



for the supermassive black holes with mass given by Eq |3.9| is 
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T = S - 



10 



,32 



-yrs 



(3.11) 



R[t,w] =/?o[l + 16Lmo(-^)f 



32LRo(-)t 
c 
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c3 
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i-r) 
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(3.7) 



which for i 



: 2LtmQ <C 1 recovers the Hulse-Taylor correc- 
tions to orbital changes given by the last term, that are linear 
in time t. In the limit of low mass stars, the dark energy cor- 
rection terms in Eqn. 13.71 are subdominant. 



R[tM = Ro 



1 MqI 
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(3.8) 



while the frequency of the emitted GW's, f — and the 
amplitude h in terms of these parameters a, P become 



/o = 



10^ 



(2P)V2, 



(3.12) 



and 



h^-lalO' 
r P 



(3.13) 



LIGO is designed to detect signals in the range / 
(100,1000)//z and amplitudes h around 10^23 
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with its optimum detection at frequencies / w 150Hz ifTsll . 
while LISA will be able to see as far as / = 10^^ Hz- One such 
example of an optimum signal for LIGO ifisll would be a black 
hole binary with parameters: a = 10, P = 10^/^,r = \25Mpc. 
Although most of the SBH's binaries, in which the dark en- 
ergy correction term is dominant, do not initially fall within 
the frequency detection limits of LIGO, they can still be de- 
tected during the time close to merging, since their frequency 
evolves as 



n = a 



(3.14) 



If, for example, we estimate the orbital changes during a 
time t ~ Hq^ with Ho the Hubble constant due to both correc- 
tion terms, then the orbital separation of a binary with initial 
frequency £2o = \Q-^^Hz, Ro = 5000pc,a = lO'^ will lately 
be n = lO^^Hz, i.e. it will fall within the detection limits. 

It can be shown that the first correction term in Eq l3.7l is 
roughly times the second correction term, therefore small 
enough to the second correction term that it can safely be ne- 
glected. 

The interesting part in Eq l3.7l is the comparison between 
the second correction term y2 which is the new term derived 
in this letter and is due to the dark energy 'bath' ,and the third 
correction term yi that has already been known and is solely 
due to the power losses from GW emission. Hereon, we in- 
clude only the second term, y2, referred to as the dark energy 
term, and the third term, yi, referred to as the GW (Hulse- 
Taylor) term. 

Replacing the values for the approriate factors in Eq |3.7| we 
obtain 



yi = 



10' 



(2P)4a 



y2 = -10"^*^(1 + w)(2p)flf ~ ±10"^''2Pflf 



39o 



(3.15) 



(3.16) 



where in the last step, the dark energy equation of state is 
taken to be (1 + w) ~ ±0.1 with the sign flipped for phantom 
dark energy, (l+w) < —0. 

For comparison, the ratio of the two correction terms is 



1045 



y2"(2P)5(i- 



(3.17) 



From Eqs l3.15H3.17l we can see that the correction term 
due to dark energy y2 can be as large as the previously known 
GW (Hulse-Taylor) correction term yi, or even dominant in 
some cases, for supermassive black holes a ^ 1 and for large 
orbital separation p ^ 2. 

Making use of this parametrization, Eqns. I3.15H3.171 we 
arrive at an important conclusion: all binaries that satisfy the 
condition 



would have merged if dark energy lies above the cc bound- 
ary, (l+w) > 0; or split apart for phantom dark energy 
(1 +w) < 0. From those, from Eqn. 13.171 all binaries with 
separation 2p > 10^ which corresponds to an initial separa- 
tion Rq > SQQOpc are dominated by the dark energy correction 
term, y2. So, Eqn. \3.18\ seems to provide a cutofffor SBH's bi- 
naries, which heavily depends on the type of dark energy and 
whether dark energy lies below or above the cc-boundary of 
w = —I. In general, measurements of /?[f,M'] can reveal the 
equation of state of dark energy (1 + w), via the Eqns. 13.161 
13.171 and binaries that fall under the category of Eqn. 13.181 
can immediately reveal whether dark energy is a phantom or 
not. For the other SBH's binaries, this information is deduced 
from 
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R[t,w] 
Ro 



[yi +y2] 



(3.19) 



2Pfl > 10** 



(3.18) 



Considering that the total mass in the universe is 10^^ kg, then 
the rough number of SBH's is around 10^ or one SBH for 
IQ-^ Mpc. If we assume that most of them are in binaries, 
then Eqn. I3.18l implies that if dark energy is a phantom then 
we should find twice as many SBH's with separation larger 
than 5000/7C as compared to the number of SBH's we would 
find if dark energy is (l+w) > for which case many of them 
would have merged to create even heavier SBH's. 

Let us take some specific examples to illustrate this effect: 
For a binary with 2p = 10^, a = 10^ that correspond to an 
orbit Ro — 10^^ m w 0.1 pc and period of rotation T = l/f = 
lO^yrs, the GW correction term is of the order yi = lO^^^f 
while the dark energy correction term is y2 = 10^^' t. In this 
example the correction term due to dark energy is many or- 
ders of magnitude larger than the GW correction and results 
in a net orbital change of 5/? ~ 10^ during a Hubble time. It 
is important to point out in this example that if dark energy 
is a phantom then y2 has a positive sign, while yi < 0. Since 
yi^yi then the merging of the two stars in the binary would 
not occur, due to the effect of the phantom energy dominat- 
ing over and compensating for the GW power losses. In fact, 
since /?[r,w] increases in this case Eq |3.7| the stars in the bi- 
nary would be 'ripped apart' over time. The lifetime of the 
stars T is equal to the Big Rip time in the phantom energy 
case and therefore larger then the characteristic times of the 
system, T and, it is larger than t^,, where is defined as the 
time when 7?[f*,w] > Rq. But, if dark energy is not a phan- 
tom or a cc A, then the dark energy corrections accelerate the 
merging process and for 2Pfl > lO'** dominate it, since in this 
case y2 > yi- By measuring the change in the orbital sepa- 
ration, using GW observations with LIGO and LISA or radio 
measurements with VLBA, we can deduce whether (l+w) is 
positive or negative, since for this class of SBH's binaries, 
the dark energy corection term y2 > yi dominates the orbital 
changes. 

However, the frequency /o = 10^"//z at initial times f = 
of this binary is not within the LIGO or LISA limits of de- 
tection. Yet, we can use X — ray or radio frequency measure- 
ments with Chandra |19] to detect changes in the phases of 
signals emitted from the SBH's binary. 
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FIG. 1: Orbital separation of the SBH's binary i^[f,w] as a function 
of time and of the dark energy's w. The parameters are chosen to be 
those of the Galaxy 0402+379 binary, with initial separation Rq ~ 
lO'^m. The range of w shown in the plot is taken: —0. 14 < [1 + w] < 
0.14. This figure shows the evolution of the binary's separation with 
time, for the case when the new corrections due to dark energy y2 
derived here, are included. 



Let us look at a more extreme class of binaries that satisfy: 
2pa « 10^**. We have Rq ~ SOOOpc, ^ 10-20(1 + w)t,yi ~ 
10-29/(2p)3f,2 ~ ±10-8 for (1 +w) ~ ±0.1. The initial fre- 
quency IS Ho ~ IQ-^^Hz and for non-phantom dark energy 
(1 + w) > 0, during f = H^^ ~ IQ-^^s we have R ~ 10"^%, 
therefore Q. ~ lO^^Hz falls within current detection limit. If 
dark energy is a phantom the orbit would be ripped apart by 
that amount, then the frequency would be extremely small and 
with no chance of detection. 

SBH's binaries have akeady been observed. We now de- 
scribe two of them which have been observed in the last 2 
years: 

i) The first example is:GALAXY 0402+379 observed in 
2007 with VLBA which has the following parameters, Rq = 
10"m,r~ 10"*i,2p== 10^fl==210^r= lO^V This binary 
thus has a frequency i2o ~ too small for LIGO de- 

tection and amplitude h ~ lO^^^. The orbital corrections due 
to GW radiation yi and dark energy accretion y2 differ by two 
orders of magnitude, yi ~ — lO^^^'^f resulting in an or- 
bital change of 5/? ~ (1 + w)10^m during each period of ro- 
tation, t = T. Since the time the signal was emitted from the 
binary, ^ j. — lO'^i, the orbit has changed by 10^'^m due 
to the GW (Hulse-Taylor) term yi, and by lO^m due to the 
dark energy accretion term y2- Thus the merging time for this 
binary, if the dark energy effect is ignored, is about another 
60,0003'r,s' since the time the emitted signal te was received. 
Included the effect of dark energy accretion results in a merg- 
ing time 1 , QOOyrs for ( 1 + w) > 0, which is roughly two orders 
of magnitude less then the standard estimate. If (1 +w) < 
then merging can not occur since y2 splitting the stars apart 
for this case, is two orders of magnitude larger than yi which 
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FIG. 2: Orbital separation of the SBH's binary -R[r,w] as a function 
of time and of the dark energy's w. The parameters are chosen to be 
those of the Galaxy 0402+379 binary, with initial separation Rq ~ 
lO'^m. The range of w shown in the plot is taken: —0. 14 < [1 + w] < 
0.14. This panel plots the same binary as the one in Fig.l but without 
the corrections induced by dark energy. The difference between the 
two plots is clear: for (1 + w) > the merging occurs faster in the 
"top" plot due to effects of y2 ; for the phantom case ( 1 + w) < 0, the 
corrections y2 slow down and even reverse the merging that the GW 
effects of the y\ term are trying to induce. 



is trying to induce merging. The evolution lifetime of each 
star in this binary is many orders of magnitude larger than the 
characteristic merging scale, T ~ ^^yrs ~ lO^^yri. The latter 
illustrates the fact that while mass changes due to dark energy 
accretion by individual SBH's are too small to give rise to de- 
tection, a binary of SBH's can yield information on dark en- 
ergy which is easily accessible by our current and near future 
experiments. 

It is quite amazing that we can obtain such a wealth of in- 
formation on the mysterious dark energy in the universe by us- 
ing local objects such as SBH's binaries, and experiments that 
were not initially designed for dark energy detection, such as 
LIGO, LISA, VLBA and Chandra experiments. This binary's 
separation as a function of time and w is plotted in Fig.l, for 
the cases when dark energy effects y2 are taken into account, 
(Fig.l); when only GW effects to the orbit, yi, are accounted 
for, but the new dark energy effects derived here are not in- 
cluded in calculating R[t\, (Fig. 2). 

ii) The second example is: a Radio Galaxy OJ287 observed 
in 2008 with VLBA |20] which is suspected to be a binary of 
two SBH's. This example is a bit trickier because one star 
is a lot heavier than the other and their total mass is 
However, although algebraically messy, it is straightforward 
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to estimate the effect of dark energy for this binary by us- 
ing the expressions of Sec. 2 for the case nii ^ 1112 and by re- 
placing Eqn. 13.21 in order to derive Eqn. 13.61 The orbital 
separation is Rq ~ or equivalently 2p ~ 10^ at f = 0. 

The current period is T ~ 123^5, and their distance from us 
r ~ 3.5Mly ~ lO^^m. If we were to approximate this bi- 
nary with equal mass SBH's each with a ~ 10^ then we get 
/o~ 10-''//z,/i~ 10-20. We receive the emitted signal from 
the distance r at a time te ~ 10 during which the orbit has 
changed by a factor {Royi) — Rol0^^^m,{RQy2) = RolO^'^m 
due to Hulse-Taylor GW effect and the dark energy effect re- 
spectively. If dark energy is ignore (y2 — 0), then merging oc- 
curs in about = lO'^iec ~ lOjOOOyri. Including the effect 
of dark energy accelerates the merging time by three orders of 
magnitude for the case (1 + w) > or stops merging and rips 
the stars apart for the case of phantom energy (1 + w) < 0. 
As mentioned, binaries are observed with the LIGO, LISA, 



VLBA, SDSS, Chandra experiments. It is possible that more 
SBH's binaries will be found in the near future by these exper- 
iments. Among other things, the SBH's binary measurements 
will allow us to distinguish more accurately whether dark en- 
ergy lies above or below the cc boundary w — ~1 and to pin 
down the value of w[z], by making use of the simple method 
proposed in this letter However, the examples discussed here 
are sufficient to illustrate the power of our method for using 
SBH binaries to obtain information about dark energy and its 
equation of state w[z] , even with existing GW and binary data, 
while avoiding background noise issues, inheritant of the large 
scale structure in the universe. 
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